Theory of dopants and defects in Co-doped Ti02 anatase 
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We report first-principles microscopic calculations of the formation energy, electrical activity, and 
magnetic moment of Co dopants and a variety of native defects in Ti02 anatase. Using these results, 
we use equilibrium thermodynamics to predict the resulting carrier concentration, the average mag- 
netic moment per Co, and the dominant oxidation state of Co. The predicted values are in good 
agreement with experiment under the assumption of O-poor growth conditions. In this regime, 
a substantial fraction of Co dopants occupy interstitial sites as donors. The incomplete compen- 
sation of these donors by substitutional Co acceptors then leads to n-type behavior, as observed 
experimentally. 

PACS numbers: PACS numbers:75.50.Pp,71.55.-i,61.72.Bb 
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I. INTRODUCTION 

The recent discovery of room temperature ferromag- 
netism in Co-doped T1O2 anastasei has led to a great 
deal of activity both to understand the origins of ferro- 
magnetic order in this material and to raise the magnetic 
ordering temperaturei^i&liSiSiLSiii One promising line of 
inquiry is to understand first the role of dopants and na- 
tive defects in the TiC>2 host. For example, the location of 
Co dopants in this material (substitutional versus inter- 
stitial), their oxidation state, and their magnetic prop- 
erties have been the subject of intense scrutinyiSi 3 * 4 *^ 
Furthermore, it has been suggested that O vacancies, 
which are believed to give rise to the observed n-type 
behavior in pure Ti02 anatase^fi may provide free elec- 
trons which mediate the exchange interaction between 
the Co dopants. 3,4 This possibility is very different from 
the hole mediated exchange interactions which are be- 
lieved to describe ferromagnetic order in a wide variety 
of other dilute magnetic semiconductors, including InM- 
nAs, GaMnAs, and MnGefii* 1 ^*^ and may have a direct 
bearing on the origin of the anomalously high Curie tem- 
perature observed in Co-doped Ti02 anatase. 

In this paper, we first use density- functional theory 
to determine the electronic structure, formation energy, 
and electrical activity of Co dopants and several native 
defects in Ti02 anatase. (In the remainder of this paper 
we will use "defects" to refer to both Co dopants and 
native defects.) We then use standard methods to cal- 
culate, as a function of temperature, the concentrations 
of each defect in the Ti02 host over a range of Co and 
O chemical potentials relevant to different growth condi- 
tions. In summary, we find that O vacancies do not play 
any significant role in Co-doped anatase. Moreover, we 
find that the observed n-type behavior in Co-doped Ti02 
strongly suggests that roughly half of the total Co con- 
tent is in interstitial sites. Under these conditions we find 
an enhancement — relative to the Co 2+ low spin state — of 
the average value of the local magnetic moment, in good 
agreement with experiment. Finally, under these condi- 
tions essentially all of the of Co — both interstitial and 
substitutional — occurs in oxidation state II, as observed 



exp er iment ally. 



II. BACKGROUND 



We begin by briefly reviewing the relevant experimen- 
tal results for Co-doped TiC"2 anatase. We concentrate 
solely on experimental results for samples in which Co 
is believed to be homogeneously distributed; hence, we 
do not address Co clustering or its consequences. We fo- 
cus on three observations which show good experimental 
reproducibility: (1) the electrical nature of the samples 
(insulating versus n-type or p-type), (2) the Co oxida- 
tion state, and (3) the average magnetic moment per Co 
dopant. 

Although the original work of Matsumoto and cowork- 
ers estimated the ferromagnetic ordering temperature to 
be larger than 400 K, the average magnetic moment 
per Co and n-type carrier concentration were modest, 
0.32 /is and 10 18 /cm 3 respectively. More recent efforts 
have led to a larger magnetic moment per Co of 1.26 
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and to n-type carrier densities of 10 /cm 
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These increased values are likely due to improvements 
in sample quality achievable with oxygen plasma-assisted 
molecular-beam epitaxy (OPMBE)£ For example, these 
samples were well characterized to rule out Co inclu- 
sions as a source of the ferromagnetismi 4 * 1 ^ Using both 
Co 2p photoemission and Co K-shell x-ray absorption 
near-edge structure (XANES), the Co dopants in TiC>2 
anatase were shown to have a formal oxidation state of 
II. A strong correlation between the magnetic and trans- 
port properties was demonstrated: both highly doped 
and highly resistive samples are typically non-magnetic, 
consistent with a picture of carrier-mediated ferromag- 
netism competing with an antiferromagnetic superex- 
change interaction^* 1 ^ 
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III. THEORY 
A. Formalism 

To understand the role of defects in T1O2 anatase, we 
initially adopt a simple picture of isolated impurities. In 
this approach, we first determine the energy required for 
the defect to form in a given charge state. We then use a 
standard thermodynamic approach to determine the ex- 
pected concentration of such defects at a given tempera- 
ture. We assume that the defects do not interact, except 
indirectly via charge transfer between them. Thus, for 
example, we do not address the origin of the apparent 
ferromagnetic coupling between Co dopants. 

The defects we consider in this work are: interstitial 
Co (Coi n t), substitutional Co on the Ti site (Co-ri), O va- 
cancies (Vo)i interstitial Ti (Tij nt ), and defect complexes 
consisting of nearest-neighbor pairs of substitutional Co 
and O vacancies (CoxiVo). There are two crystallo- 
graphically distinct types of such complexes: a CoTiVo 
pair oriented along the c-axis and a CoxiVo pair oriented 
nearly in the ab plane. We refer to these as CotiVo-c and 
CoTiVo-a&, respectively. 

We make no assumptions about which defects are 
donors and which are acceptors, whether they are neu- 
tral or charged, or how many are actually present. At a 
given temperature, the concentration of each defect (in a 
given charge state) will be determined by its formation 
energy. We use density-functional theory in a supercell 
approach to calculate these formation energies according 
to 

EL m = Ef-J2nm+QE F , (1) 

j 

where is the total energy of a supercell containing one 
defect in charge state q; rij and /ij are the number and 
chemical potential (energy per atom in the ground-state 
phase) of each atomic species in the supercell; and Ep is 
the Fermi energy, measured with respect to the valence 
band maximum (VBM) of the hosti 16 i 17 i 18 i 19 i 20 Although 
the Fermi energy appears to be an independent variable 
here, it is in fact determined by the constraint of elec- 
troneutrality, as described below. Also, we note that al- 
though the total energies on the right-hand side of Eq. 2] 
depend on whether all-electron or pseudopotential meth- 
ods are used, the formation energy itself is analogous to 
a binding energy, and can therefore be accurately calcu- 
lated within either method. 

In thermal equilibrium, the concentration of each de- 
fect, D, is determined by its formation energy: 

C q D =N sites exp(-El rm /k B T), (2) 

where N s i tcs is the number of sites per unit volume avail- 
able to the defect. Since the formation energies depend 
on the chemical potentials and the Fermi level, it is evi- 
dent that the concentrations also depend on these quanti- 
ties. The concentrations must also satisfy the constraint 



of overall electrical neutrality; this provides an additional 
equation that we use to determine Ep for any given 
choice of chemical potentials. The electroneutrality con- 
dition must take into account the contributions not only 
from charged defects, but also fromp and n (the hole and 
electron carrier densities). Thus, for each given choice of 
the oxygen chemical potential, [Ao, and cobalt chemical 
potential, fico, the following equation must be numeri- 
cally solved: 

p(E F ) - n(E F ) +Y,qC q D (E F ;fio^Co) = 0. (3) 

D,q 

Here, the sum is over all defects, D, in all possible charge 
states, q. The Ti chemical potential does not enter ex- 
plicitly into this equation for reasons discussed in the 
following section, and, of course, the Co chemical poten- 
tial is only relevant for defects involving Co. The carrier 
densities p and n are evaluated using the conventional 
semiconductor expressions along with the ab initio den- 
sity of states of the host material, with a scissors operator 
applied to give the experimental band gap. 

B. Chemical Potentials 

The atomic chemical potentials, Hj, on the right hand 
side of Eq. ^ are closely related to the experimental 
growth conditions. A high value of chemical potential 
of a particular atomic species is equivalent to a growth 
environment that is rich in that species (in the sense of 
high partial pressure). 

The chemical potentials of Ti and O which enter into 
Eq. n are not independent: equilibrium between the Ti 
and O atomic reservoirs and bulk Ti02 anatase requires 
that 

Mti + Mo 2 = MTio 2 , (4) 

where /XTio 2 is the total energy of bulk anatase. More- 
over, in order to preclude the precipitation of bulk hep 
Ti or O2 dimers there are additional thermodynamic re- 
strictions on the individual chemical potentials: 

Mo < Mo; ( 5 ) 

and 

MTi < M^' k , (6) 

where /Xq is the energy per atom in an O2 dimer in its 
spin triplet ground state and /i^" 11 * is the energy per atom 
in hep Ti. It is conventional to refer to the upper limit 
in Eq. [31 as the O-rich limit and, similarly, to the upper 
limit in Eq. [5] as the Ti-rich limit. Because of the rela- 
tionship between and /io in Eq. there are only two 
independent variables in this approach: the O and Co 
chemical potentials. This explains why the Ti chemical 
potential does not enter explicitly into Eq. [3J 
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Thus, for a given choice of /^co and /^o, all the other 
quantities (including Ep) in Eq. [3|are completely deter- 
mined. In practice, we eliminate /ico as an independent 
variable by constraining the total Co concentration to 
a typical experimental value (5%). Hence, in our for- 
mulation, the concentration of all defects is entirely de- 
termined by the choice of O chemical potential and the 
constraint of total Co concentration. 



C. Temperature 

Although temperature appears in both Eq. ^ (implic- 
itly) and Eq.0 (explicitly), our results are not especially 
sensitive to this variable. Moreover, we stress that the 
temperature appearing in Eq.|21has no connection to the 
magnetic ordering temperature, and should instead be 
understood simply as the temperature at which we evalu- 
ate the concentration of defects in equilibrium with their 
respective reservoirs. Upon completion of growth, we 
consider the reservoirs to be disconnected and hence the 
number of these constituents to be fixed. Thus we set 
the temperature to a particular value and examine the 
behavior of the system as a function of the O chemical 
potential. To this end, we calculate the defect concentra- 
tions using a typical growth temperature of 873 Ki2i£ 



D. Computational Details 

The total energy calculations were performed in a su- 
percell consisting of a 3x3x2 periodic repetition of the 
primitive unit cell; thus, for the pure TiC>2 host, these 
supercells contain 108 atoms. Structural relaxation was 
performed within a sphere of radius of 4 A centered on 
the defect in question; calculations using a sphere radius 
of 5 A give the same total energy to within 50 meV. This 
approach has proven both efficient and accurate for other 
similar metal-oxide insulators ^* The total energy calcu- 
lations were evaluated in the local-density approximation 
(LDA) within the ultrasoft pseudopotential formalism 21 
as implemented in the VASP code^ with the zone cen- 
ter used to sample the Brillouin zone of the supercells. A 
kinetic-energy cutoff of 400 eV was used in all total en- 
ergy evaluations. Only defects involving Co atoms were 
treated in a spin-polarized fashion. The ground state of 
bulk hep Ti was treated non-magnetically, bulk hep Co 
was treated in an ferromagnetic configuration, and the 
total energy of the 02 reference dimer was evaluated for 
the spin triplet. 

Since the Fermi energy in Eq. ^ is measured with re- 
spect to the VBM of the host, we must align the VBM 
in the charged defect supercell with that of the host ma- 
terial. To this end, we used a local site average of the 
electrostatic potential to define a reference energy; this 
local average was evaluated by a test charge approach 
in which we calculate the electrostatic energy of a nar- 
row Gaussian charge distribution far from the defect. We 



have checked that using, as an alternative, the Ti 3p semi- 
core eigenvalues gives very similar results. Including the 
zone sampling, kinetic-energy cutoff, lattice relaxation 
and choice of reference energy, we estimate the numerical 
uncertainty in our results to be 100-150 meV, sufficient 
for addressing trends with respect to growth conditions. 

To compute the total energy of a charged periodic sys- 
tem, we use a standard procedure which systematically 
corrects for the artificial and slowly converging Coulomb 
interaction between charged defectsi2&24 A neutralizing 
homogeneous background charge density is first added 
to the total charge density; this makes the total energy 
a well-defined quantity. Next, the artificial interaction 
among the charged defects within this neutralizing back- 
ground is subtracted from the total energy; this interac- 
tion is estimated by expanding the defect-induced elec- 
tron density in a multipolc series up to quadrupole order 
and then computing the contribution to the total energy 
from the interaction of these multipoles^ 

Finally, we note that although the electronic band gap 
does not enter explicitly into Eq. the LDA underes- 
timation of the gap does affect the formation energies 
of shallow donors, such as O vacancies and interstitial 
Co. We have observed, however, that the depths of these 
donor levels (relative to the conduction-band edge) are 
essentially independent of the value of the band gap. 
We demonstrated this by calculating the position of the 
donor level using an artificially reduced lattice constant; 
as the band gap increases with decreasing lattice con- 
stant, the donor level closely tracks the conduction-band 
edge. Thus, we have corrected the formation energies of 
these donors using the experimental value of the band 
gap and the depth of the donor levels given within LDA. 
In addition, we assume that the LDA correctly predicts 
the formation energies of the +2 charge states of these 
defects, because all defect levels are empty in this charge 
state. 



IV. RESULTS AND DISCUSSION 
A. Electronic Structure of Coti 

We discuss first the electronic structure of the isolated 
substitutional Co dopant. Fig. Qfa) shows the single- 
particle levels of neutral substitutional Co, as determined 
from the eigenvalue spectra at k = T using a 108-atom 
supercell. The point group for this defect is D2d, so that 
all but one of the 3<i-orbital degeneracies are lifted by the 
crystal field; nevertheless, the environment surrounding 
the Co dopant is still very nearly cubic, so that the e g 
and t2 g parentage of these levels is easily seen. There are 
3 occupied majority-spin levels and 2 occupied minority- 
spin levels, yielding a net magnetic moment of 1 \xb- 

It is instructive to study separately the exchange split- 
ting and crystal-field splitting. This separation occurs 
naturally in the q = — 1 charge state, for which the ex- 
change splitting vanishes, as shown in Fig.^b). In this 
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case the largest crystal-field splitting, about 1 eV, is be- 
tween the e g and t 2g manifolds, with the e g manifold 
near the LDA conduction-band edge and the t2 g near the 
valence-band edge. The Did crystal field further splits 
the e g manifold into an upper a\ level of z 2 symmetry 
and a lower b\ level of x 2 — y 2 symmetry; this splitting is 
0.3 eV. Likewise, the crystal field splits the ti g manifold 
into an upper hi level of xy symmetry and a two-fold e 
level of (xz, yz) symmetry; this splitting is 0.2 eV. 

For charge states with non-zero magnetic moments, we 
find that the exchange splitting varies approximately lin- 
early with the magnetic moment. (Since the orbital mo- 
ment is strongly quenchedj2& we regard the magnetic and 
spin moments as equivalent.) For neutral substitutional 
Co, the exchange splitting within the e g manifold is 0.2 
eV for both levels. Within the t2 g manifold, the exchange 
splitting is strongly orbital-dependent: 0.5 eV for the 62 
levels and 0.3 eV for the e levels. 

The calculated magnetic moments of substitutional Co 
in different charge states, summarized in Tabled can now 
be easily understood from the results of Fig. Q For the 
neutral substitutional, there is a two-fold degenerate half- 
filled minority spin level at the Fermi level. Removing an 
electron from this level leads to a +1 charge state with 
a moment of 2 \ib- Likewise, adding an electron to this 
level leads to a —1 charge state with zero moment, since 
the ti g manifold is now completely filled. 

The stable charge states and magnetic moments for 
interstitial Co are quite different from the substitutional 
case. For example, we find that all charge states experi- 
ence off-center structural relaxations of order 1 A, lifting 
all remaining orbital degeneracies. The +2 charge state 
of interstitial Co, which is the lowest-energy charge state 
over most of the gap, has a moment of 1 \lb- The +1 
charge state leads to a moment of 2 fj,s, whereas the 
neutral configuration has a moment of 1 [Ib- 



B. Oxidation State 

Since the correspondence between "formal oxidation 
state" and "charge state" will be an important issue here, 
we briefly summarize the relationship between the two. 
A similar discussion has been given for transition-metal 
impurities in GaP. 1 - When a neutral Co substitutes for 
Ti in Ti02, it takes on the oxidation state IV, the same as 
that of Ti. Hence a substitutional Coxi in the —1 charge 
state has a formal oxidation state of III, and a substi- 
tutional Coxi with charge —2 has an oxidation state of 
II. In a similar fashion one can relate the oxidation state 
and charge state of interstitial Co; in this case however, 
the charge state of the dopant is the same as the oxida- 
tion state. Hence, neutral interstitial Co has oxidation 
state 0, the charge state +1 has oxidation state I, and so 
forth. The oxidation state of Co-related dopants in their 
various stable charge states are summarized in Table [I] 
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FIG. 1: First-principles energy-level diagram for substitu- 
tional CoTi in Ti02 anatase. (a) Neutral charge state; (b) 
Singly charged state, q = — 1. The density of states of the 
anatase host is shown as the gray shaded regions with a scis- 
sors operator applied to give the experimental band gap. The 
length of each impurity level is proportional to its d character. 
The exchange splitting, Ax, and crystal field splitting, Acf, 
are discussed in Sec. 4(c) (see text). 



TABLE I: Stable charge states of various Co dopants in Ti02 
anatase, with their formal oxidation states and calculated 
magnetic moments. The results for the —2 charge state of 
CoTi are obtained from a model described in the text. The 
oxidation state values in parentheses for CoTiVo complexes 
are based on the assumption that the neutral complex can be 
represented as CoTi and Vo in the —2 and +2 charge states, 
respectively. 
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C. The 2 Charge State of Co Ti 

In the LDA calculation there is no stable —2 charge 
state of substitutional Co: upon adding an additional 
electron to the stable —1 charge state, we find no dopant- 
derived level within the gap, but rather partial occupa- 
tion of the LDA conduction bands. This is problematic, 
since the experimental finding of an oxidation state of 
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II for Co would be consistent with a —2 charge state 
for substitutional Co. Here we investigate whether the 
absence of a —2 charge state within LDA is due to the 
well-known underestimation of the band gap by the LDA 
(the LDA predicts a band gap of 2.2 eV compared to the 
experimental value of 3.2 e\Sl). 

Rather than attempting to correct the LDA band gap, 
we propose a more direct description and try to estimate 
the formation energy of the —2 charge state using our 
results for the energy levels of the —1 charge state. Since 
the formation energy is in general a linear function of 
the Fermi energy, we do this in terms of the "energy of 

transition," E(—J ), which is defined as that value 

of E F for which E^~~ is equal to Ef~~ . Referring to 
Fig. W[b), we model the —2 charge state by occupying 
the first available impurity level in the q = — 1 level di- 
agram, namely the empty b\ level. In this diagram, the 
bi level is at Ep + Acf, so that occupying it will shift 
the Fermi level upward by an amount Acf ■ However, by 
singly occupying this level, one expects an accompany- 
ing exchange splitting and thus a downward shift of the 
Fermi level. Hence, we estimate the change in the Fermi 
energy between the —1 and —2 charge state to be given 
by 




E F (eV) E F (eV) E p (eV) E p (eV) 

FIG. 2: Formation energies of (a) Coxi, (b) Coi nt , (c) 
CotiVo, (d) Vo defects as a function of Fermi level in the 
O-rich limit. The Co chemical potential was chosen to give a 
total Co concentration of 5%. The Fermi energy which sat- 
isfies electroneutrality for these choice of growth conditions 
and Co concentration is denoted by the dotted vertical line. 
The stable charge states for each defect are labeled. 



E(-/--) = E(0/-) + A CF -A x , (7) 

where E(0/— ) is the energy of transition between the 
neutral and —1 charge states; Acf is the crystal- field 
splitting between the b 2 and b\ levels; and Ax is the 
exchange splitting due to the unpaired b\ electron. 

We approximate Ax by the b\ exchange splitting for 
the neutral state; this is 0.3 eV, as shown in Fig^a). We 
neglect the onsite energy which arises from the interac- 
tion of the added b\ electron with the ti g electrons, but 
we expect this energy to be significantly smaller than the 
crystal- field energy, and hence expect Eq.[7|to be reason- 
ably accurate. 

A similar proposal could be made for more-negative 
charge states, for example, the —3 charge state of sub- 
stitutional Co. However, we should then include an ad- 
ditional term, U , for the onsite interaction between two 
b\ electrons of opposite spin. This onsite energy will be 
large and can be estimated as 

U > E (0/-) - E(+/0) - A x , (8) 

where Ax = 0.3 eV is again taken from the neutral con- 
figuration, but in this case accounts for the exchange 
energy lost when the moment is reduced from 1 /is 
to zero in the transition from the neutral to the — 1 
charge stated We find U > 1.0 eV and, hence, the 
—3 charge state of substitutional Co will lie well above 
the conduction-band edge, so that it (and more highly 
charged negative states) can be excluded from further 
consideration. 



D. Formation Energies in the O-rich Limit 

To address the characteristics of Co-doped samples 
grown with OPMBE we first consider conditions O-rich 
growth conditions. Figure [21 shows the defect formation 
energies versus Fermi energy in the O-rich limit, with the 
total Co concentration constrained to be 5%. In compari- 
son to experiment this scenario has three main shortcom- 
ings: (1) the average magnetic moment is 1 /is, signifi- 
cantly smaller than the measured valued of 1.26 [Ib] (2) 
the Fermi level is well below midgap, whereas experimen- 
tal measurements show the material to be K-type^^ 
(3) Co appears almost entirely as neutral substitutionals 
(oxidation state IV), whereas Co 2p photoemission and 
XANES suggest the oxidation state of Co is IlAii These 
discrepancies suggest that the conditions present during 
growth of the samples are not O-rich. 



E. Variation with O Chemical Potential 

To determine what growth conditions give rise to the 
observed magnetic and transport properties of Co-doped 
anatase samples, we study the consequences of varying 
the O chemical potential away from its upper limit. In 
Fig. |3 we plot, versus the O chemical potential, the fol- 
lowing quantities: concentration of of substitutional and 
interstitial Co; n and p carrier densities; average mag- 
netic moment per Co; percentage of Co in different ox- 
idation states; and concentration of O vacancy-related 
defects (the concentration of interstitial Ti is negligible 
throughout this range of O chemical potential). For every 
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FIG. 3: Variation, versus O chemical potential, of (a) CoTi 
concentration, (b) Coint concentration, (c) carrier concentra- 
tion, (d) average magnetic moment per Co, (e) percent of 
Co in various oxidation states, and (f) concentration of Vo- 
related defects versus O chemical potential. In panels (a) 
and (b) the concentration of each charge state of Co is indi- 
cated as well as the total concentration. In panel (f) only the 
concentration of CoTiVo defects are shown, since the concen- 
tration of isolated O vacancies is much smaller. Values of the 
O chemical potential to the left of the dashed vertical line 
are those for which our theoretical results agree qualitatively 
with experiment (see text). 



value of the O chemical potential, the total concentration 
of Co was constrained to 5%. 



1. Cobalt concentration 

For values of the O chemical potential within ~1.5 eV 
of the O-rich limit, the substitutional site remains the 



preferred Co site. The situation is very different in the O- 
poor limit. In this regime, interstitial Co (in the +1 and 
+2 charge states) becomes energetically competitive with 
substitutional Co, and hence the concentrations of sub- 
stitutional and interstitial Co become comparable. This 
results from a subtle balance of the O and Co chemical 
potentials. As one decreases the O chemical potential 
from its upper limit, the formation energies of substi- 
tutional dopants increases, since it costs more energy to 
replace Ti with Co as one moves toward the Ti-rich limit. 
Of course, as a function of the O chemical potential, 
the Co chemical potential must also change, increasing 
monotonically to maintain the fixed total Co concentra- 
tion. At ~1 eV below the O-rich limit, interstitial Co 
begins to play a role; the preference for the interstitial 
site increases monotonically below this point, eventually 
accounting for most of the total Co concentration in the 
O-poor limit. 



2. Carrier densities 

The defect concentrations shown in Fig. [3] were com- 
puted for a temperature of 873 K, as mentioned earlier. 
However, since transport, magnetometry, and photoemis- 
sion measurements are generally performed at room tem- 
perature, we have used a more relevant temperature of 
300 K to compute the carrier concentrations, magnetic 
moments, and oxidation-state fractions — while keeping 
the total defect concentrations themselves frozen at their 
high-temperature values. 

Figure [^Jc) shows the logarithm of the hole and elec- 
tron carrier densities evaluated at 300 K as a function of 
the O chemical potential. Near the O-rich limit the car- 
rier densities are nearly constant, since the Fermi level 
maintains a value of ~0.8 eV relative to the VBM due 
to compensation by equal but small numbers of substitu- 
tional Co in the —1 and +1 charge states. In this region 
the material is very weakly p-type with a hole concentra- 
tion ^10 6 /cm 3 and will likely appear insulating in trans- 
port measurements. 

At ~1 eV below the O-rich limit, the appearance of in- 
terstitial Co leads to partial compensation of the Co sub- 
stitutional defects. This compensation drives the Fermi 
level towards the conduction-band edge, and leads to a 
marked increase in the electron concentration as the O 
chemical potential is further reduced. The step-wise be- 
havior of the carrier density as the O chemical potential 
is decreased results from the Fermi level passing through 
donor levels due to interstitial Co. At the O-poor limit, 
we find n ~ 10 20 /cm 3 . Thus, we suggest that the ex- 
perimentally observed n-type behavior results from the 
incomplete compensation of interstitial Co by substitu- 
tional Co, which facilitates the thermal excitation of elec- 
trons from interstitial Co into the conduction band. 
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3. Magnetic moment 

In Fig.|3Jd) we show the average magnetic moment per 
Co, defined here as 

D,q D' ,q' 

where M q D is the magnetic moment of a Co-related defect 
in charge state q. We assume a ferromagnetic alignment 
of all the moments, and thus M can be considered an up- 
per bound for the measured average value of the moment 
per Co. 

The average moment shows a non-monotonic variation 
with O chemical potential. In the O-rich region it is con- 
stant with a value of 1 hb, because neutral substitutional 
Co is the dominant defect. For intermediate values of the 
O chemical potential, M is less than 1 due to the pres- 
ence of substitutionals in the —1 charge state, which have 
zero moment. In the O-poor limit the average moment 
per Co is larger than 1 /x_b, due to the appearance of in- 
terstitial Co in the +1 charge state, which has a moment 
of 2 [Ib ■ Hence, only in the O-poor limit do we obtain a 
moment per Co larger than the Co 2+ low-spin value of 1 
fiB, and thus consistent with experiment. 

4- Oxidation states 

Figure HJe) shows the percentage of Co dopants in dif- 
ferent oxidation states. In the O-rich limit the oxidation 
state of Co is IV, since it occurs primarily as a neutral 
substitutional. This oxidation state dominates until ~1.5 
eV below the O-rich limit where oxidation state III, re- 
sulting from substitutionals in the —1 charge state, is 
briefly dominant. Between 3 and 4 eV below the O-rich 
limit the predominant oxidation state of Co is II, the 
same as that deduced from the 2p photoemission and 
XANES results of RefsQ andQ In the O-poor limit the 
oxidation state is primarily I, since the dominant type of 
defect is interstitial Co in the +1 charge state. Although 
this oxidation state has not been observed experimen- 
tally, comparison to photoemission or XANES results on 
samples with known Co oxidation state of I has not been 
established. 



5. O vacancies 

For all thermodynamically allowed chemical potentials, 
the concentration of isolated O vacancies is negligible, as 
well as that of Co-riVo complexes; their concentrations 
are shown in Fig. |3Jf ) to be never above 10~ 5 %. These 
low concentrations result from the fact that the forma- 
tion energies of these two defects are much higher than 
the thermal energy fc^T = 60 meV; this is clear from 
Fig. |21 where the formation energy of both isolated Vo 
and CotiVo complexes is greater than 2 eV for any value 



of the Fermi level. Thus O vacancies, either isolated or in 
complexes with substitutional Co, will play no significant 
role in determining the carrier concentration in Co-doped 
samples. 

We note that the formation energies of O vacancies in 
Fig. |21 and the deduced donor level positions in the band 
gap are incompatible with the interpretation of the ob- 
served n-type conductivity of pure Ti02 samples^ Even 
in the O-poor limit, the LDA does not lead to a car- 
rier density with the observed value 10 18 /cm 3 , nor does 
it give the same temperature dependence of the carrier 
concentration, as the energy for activation is quite differ- 
ent: ^200 meV (the position of the highest O vacancy 
donor level below the conduction band edge) in the LDA 
versus 4.2 meV in the results of Ref. We do not have 
a definitive resolution to this apparent disagreement; we 
have checked, however, that if we adjust the O vacancy 
formation energies to reproduce the experimentally ob- 
served carrier density and temperature dependence in 
undoped samples, the conclusions of this work are un- 
changed. 

Regarding the experimental interpretation that O va- 
cancies are also the source of n-type carriers in Co doped 
samples we note that the observed average moment 
per Co of 1.26 /xs can only be reproduced if some fraction 
of Co interstitials in the +2 charge state are present, as 
only these defects have a moment larger than 1 fis ■ Nei- 
ther substitutional nor interstitial Co has any significant 
orbital moment^ so that the observed magnetic moment 
per Co can only result from the statistical distribution of 
Co moments in interstitial (2 hb) and substitutional (1 
Hb) sites. In such a situation the Fermi level and carrier 
concentration are determined solely by the Co dopants. 



6. Conclusion 

Considering these trends in the n-type carrier density, 
average magnetic moment per Co, and oxidation state 
of Co, we suggest the actual conditions of growth of 
Co-doped Ti02 anatase are O-poor, corresponding to O 
chemical potentials between 3.8 and 4.6 eV below the O- 
rich limit (values of the O chemical potential to the left 
of the dashed vertical line in Fig. [SJ. Under these type 
of conditions we obtain qualitatively and quantitatively 
good agreement with experimental observations. 



V. SUMMARY 

In summary, we have examined the role of native 
defects and Co dopants in Ti02 anatase over a range 
of chemical potentials corresponding to different growth 
conditions. Under O-rich growth conditions we find that 
Co dopants will be formed primarily in neutral substitu- 
tional form corresponding to oxidation state IV, an aver- 
age magnetic moment of 1 [Ib and insulating electrical 
character. These results are in conflict with the exper- 
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imentally observed behavior of Co-doped samples and 
suggest that the growth conditions are more likely to be 
O-poor. O-poor conditions lead to roughly equal con- 
centrations of substitutional and interstitial Co, n-type 
behavior resulting from thermal excitation of electrons 
from interstitial Co into the conduction band, and an av- 
erage magnetic moment per Co in good agreement with 
experiment. 
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